Transverse Ising Chain under Periodic Instantaneous Quenches: Dynamical 
Many-Body Freezing and Emergence of Slow Solitary Oscillations 



Sirshendu Bhattacharyya 
R.R.R. Mahavidyalaya, Radhanagar, Hooghly, India 

Arnab Das* 

Theoretical Division (T-4), LANL, MS-B213, Los Alamos, New Mexico - 87545, USA and 
Max-Planck Institute for the Phystcs of Complex Systems, 
Nothnitzer Str. 38, Dresden 01187, Germany (present address). 

Subinay Dasgupta 

Department of Physics, University of Calcutta, 92 Acharya Prafulla Chandra Road, Kolkata 700009, India 

(Dated: August 13, 2012) 

We study the real-time dynamics of a quantum Ising chain driven periodically by instantaneous 
quenches of the transverse field (the transverse field varying as rectangular wave symmetric about 
zero). Two interesting phenomena are reported and analyzed: (i) We observe dynamical many-body 
freezing (DMF) [l], i.e. strongly non-monotonic freezing of the response (transverse magnetization) 
with respect to the driving parameters (pulse width and height) resulting from equivocal freezing 
behavior of all the many-body modes. The freezing occurs due to coherent suppression of dynamics 
of the many-body modes. For certain combination of the pulse height and period, maximal freezing 
(freezing peaks) are observed. For those parameter values, a massive collapse of the entire Floquet 
spectrum occurs, (ii) Secondly, we observe emergence of a distinct solitary oscillation with a single 
frequency, which can be much lower than the driving frequency. This slow oscillation, involving 
many high-energy modes, dominates the response remarkably in the limit of long observation time. 
We identify this slow oscillation as the unique survivor of destructive quantum interference between 
the many-body modes. The oscillation is found to decay algebraically with time to a constant value. 
All the key features are demonstrated analytically with numerical evaluations for specific results. 



I. INTRODUCTION 

Dynamics of driven quantum many-body system is 
an emerging paradigm for studying and unveiling new 
quantum phenomena. Last few years have witnessed a 
surge of theoretical endeavors in understanding dynamics 
of quantum many-body systems under simple drivings. 
A major part of these recent activities is concentrated 
around quantum quenches, leading to several interesting 
and novel issues including (but not limited to) univer- 
sal quench dynamics across quantum critical points - as- 
sociated quantum Kibble-Zurek mechanism, physics of 
non-equilibrium excitations, and the physics of thermal- 
ization in quantum systems (see for a review, Ref. |2^, ; and 
C. de Grandi et. al., S. Mondal et. al., and U. Divkaran 
et. al, in Ref. [5] and references therein). The main fo- 
cuses of these studies, e.g., the final defect density in a 
quantum quench, or the effective temperature in a ther- 
malized system, however, are insensitive to the details 
of the quantum coherence of the underlying many-body 
dynamics. For example, the dynamical idea behind quan- 
tum Kibble-Zurek mechanism (4j is a robust translation 
of the classical Kibble-Zurek idea [5] to quantum systems 
- of course, the origin of the relevant length scales and 
time scales are different. 

Here we focus on another important class of driven 
quantum non-equilibrium phenomenon, where quan- 
tum coherence plays the central role. Though non- 
adiabaticity is a common covering for all interesting non- 



equilibrium phenomenon, here coherent quantum me- 
chanical suppression of dynamics contributes crucially to 
the non-adiabaticity of the dynamics which makes the re- 
sulting response behavior difficult to explain using classi- 
cal intuitions. We discuss dynamics of periodically driven 
quantum many-body system. Coherent periodic driv- 
ing can give rise to surprising phenomenon in quantum 
many-body system, that counters our classical intuitions 
drastically [1 . The role of quantum coherence in the im- 
portant context of superfluid-insulator transition realized 
in periodically driven optical lattice was demonstrated 
earlier [6l[3. Owing to the experimental break-through 
in attaining long coherence time in quantum many-body 
systems in last decade, for example, within the frame- 
work of atoms/ions in optical lattices and traps, this 
coherent regime is becoming more and more accessible 
experimentally (see, e.g., [SHU])- Here we study the co- 
herent dynamics (Schrodinger dynamics at zero temper- 
ature) of a simple paradigmatic system - the transverse 
Ising chain [13^ subjected to a train of rectangular pulses 
of the transverse field. Two interesting phenomenon are 
reported - both purely quantum mechanical in origin and 
are results of coherent many-body dynamics. 

It has been observed recently that a class of inte- 
grable quantum many-body systems exhibit the phe- 
nomenon of dynamical many-body freezing (DMF), i.e. 
non-monotonic freezing behavior of all the many-body 
modes when driven externally by varying a parameter 
in the Hamiltonian continuously 1 . The said freezing 
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behavior is counter-intuitive to the "classical" picture of 
a driven system falling out of equilibrium. The classi- 
cal behavior arises from competition of two timescales: 
the driving period and the relaxation time of the sys- 
tem (see, however, [2]). In contrary to the expected 
monotonically increasing freezing behavior of the system 
with respect to the driving frequency according to that 
picture, we observe strongly non-monotonic freezing be- 
havior, with maximal freezing for certain combination of 
driving amplitude and frequency. A related phenomenon, 
observed in context of a single particle localized in a pe- 
riodically driven potential - known as dynamical local- 
ization, or synonymously, coherent destruction of tunnel- 
ing (CDT) is weU studied pMl9]. In [1^, interestingly, 
it has been shown in the context of periodically driven 
BEC, that the driving can lead to steady BEC-like states 
which are different from the equilibrium ground state of 
the undriven Hamiltonian. Above findings motivate us 
investigating such phenomenon in a pulse driven many- 
body system, where, instead of a smoothly varying driv- 
ing rate, we have sequences of instantaneous quenches 
and subsequent waiting times. Here we observe DMF, 
confirming the generality of the phenomenon beyond si- 
nusoidal driving. We deduce the exact condition for the 
maximal freezing analytically and explore other charac- 
teristics of the freezing phenomenon. 

In addition to DMF, we observe another interesting 
phenomenon away from the freezing peaks. In the limit of 
long observation time, we see spectacular dominance of a 
single long-lived oscillation (with frequency much smaller 
than the driving frequency) in the response dynamics. 
Surprisingly, this happens even in the limit of strong and 
fast driving (pulse amplitude and frequency much larger 
than the inter-spin coupling). We discuss the origin and 
nature of this intriguing quantum oscillation. 



II. THE MODEL AND THE DYNAMICS 

We quench the transverse field F from +Fo to — Fq 
and back in successive time intervals of duration T in a 
transverse Ising chain Hamiltonian: 



N 
.7 = 1 



N 



(1) 



where the field F(t) varies like a square- wave with period 
T at i = : 



Fo for nT <t< in+^)T 
-To for {n+^)T <t<{n+ 1)T 



(2) 



with n = 0, 1, 2, • • • and Fq > 0. We set the energy scale 
by taking J — 1. In order to investigate the dynamics in 
this case, first we diagonalize Hamiltonian ([!]) for a given 
value of F by Jordan- Wigner transformation followed by 
Fourier transform f20j . This transforms the Hamiltonian 
([l]) into a direct-sum of Hamiltonians of non-local free 



fermions of momenta k. The Hamiltonian preserves the 
parity of the fermion number (even/odd) and the ground 
state always lies in the even-fermionic sector. We work 
with the projection of the Hamiltonian in this sector, 
given by 



n 



fe>0 

= (— 2isinfc) 
2(F + cosfc) 



(3) 



aia 



t 
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where fc= (2n + l)7r/A^; n = 0, 1, A^/2- 1. The ground 
state of T-Lk is a linear combination of the fermionic oc- 
cupation number basis states |0)fe = |Ofe,0_fe) (both ±k 
levels unoccupied) and |l)fe = |lfe,l_fc) (both ±fc levels 
occupied), and the Hamiltonian does not couple them 
with the two other basis states |Ofe,l_fc) and |lfe,0_fc). 
Hence starting with ground state, the dynamics always 
remains confined within a manifold which is the direct 
product of the 2— dimensional subspaces spanned by |0)fc 
and We denote the eigenstates of T-Lk within this 

subspaces as |(F, fc)_) (ground state), and |(F, A:)+), with 
eigenvalues — A(F,fc), A(F,fc), where 



A(F,fc) = 2v/F2 + l + 2Fcosfc 
|(F,fc)_) = ^cos0|l)fc-sin0|O)fc 
|(F,fc)+) = ^sin0|l)fe+cos0|O)fc 

— sin k 

tan = 



F + cos fc + \/F2 + 1 + 2F cos /fc 
We now solve the Schrodinger equation 

d\i^k) 



ih- 



dt 



'Hk\'ipk), 



(4) 
(5) 
(6) 

(7) 



(8) 



where the wave-function in the time-dependent energy 
eigen-basis may be expressed as 

l^fe) =a;_(t)|(F,fc)_(<))+x+(t)|(F,fc)+(t)). (9) 

If T{t) is constant (say, Fq) over a time interval to i, 
then we have 

a;±(i) = x±(io)exp|T-^(t-to)A(Fo,fc)|. (10) 

At time t = let the system be in a state 

iVfe) =«|(Fo,fc)_)-f /3|(Fo,fc)+) (11) 



with \aY + = 1- Then according to Eq. (10 1 aX t ^ 



^ — e (where e is a small positive number) , the coefficients 



are given by. 



where 



e) 







T 
2h 



A(Fo,fc). 



(12) 



(13) 
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Using the continuity of l^fk) sd, t — one obtains the 
wave function at t = + e in terms of |(— Fq, fc)-) and 
|(— ro,fc)+). Time evolution in the second half proceeds 
in the same way as in the first half and the transformation 
over one full cycle is given by, 



x_{T + e)\ _ fa 
x+{T + e) -^"[P 



where. 



cos (p — sm ( 
sin (f> cos c 

cos (f) sin (/) 
— sin (h cos (b 



e~ 
e- 



(14) 



(15) 



= |^A(-ro,fc) and = 01-02 (16) 



Here 



where 9i, 02 are the values of 6 (as defined in Eq. ([7|)) 
for r = +ro, —To respectively. 

During the first half-cycle after n full cycles, at a time 
t = nT + r with < r < , the coefficients are given by, 



x^(nT + t) 
x+{nT + t) 



eV3 




(17) 



where /ia — j^X{Tq, k). Similarly, during the second half- 
cycle after n full cycles, at a time t = nT + + t, the 
coefficients are given by. 



X- (nT + ^ + T} 



T 



x+{nT+ ^ +t) 
cos (f) sin 4)\ ( e^^^ 



gV4 

e-^'"* 



(18) 



sm cos ( 







/3 



where /i4 = ^A(— Fq, k). 

Transverse magnetization (per spin) at any time is 
given by. 



i+iy-Mfc=-^i+- r 

N ^ TT Jo 



Mk dk (19) 



where Mk = ^{tpk \ {alak + a.*_f.a-k) iV'fe)- From Eqs. E 



Mfc — I (a;_ cos0j -I- Xj^ sinOj) Y 



(20) 



with j = 1,2 according as we are in the first or second 
half-cycle respectively. 

In order to calculate UjJ , giving the time-evolution af- 
ter nth cycle, we note that for any 2x2 matrix, 

U2 =-(TrU,)l + (detUfc)Ufc 

This shows that one can write 

UI! = a„l + fo„Ufc (21) 



The recursion relations for a„ and 6„ can be easily solved 
to get 

a„ = — &ri-i and 6„ — sm(nu)k)/ sinw^. (22) 

where coswfc = cos(/^i -f ^2) cos^ -I- cos(/^i — ^2) sin^ 4>. 
The expressions 5„ are the Chebyshev polynomials of the 
second kind in coscj^. 



III. DYNAMICAL MANY-BODY FREEZING 
(DMF) 

The system is initially {t — 0) in the ground state of the 
Hamiltonian with F = +Fo, before it is driven by the 
pulses. We have computed the magnetization numeri- 
cally at any time (within a cycle) by obtaining x_ and 
x+ from Eqs (ITtI [T9| , substituting them in Eq. (pO| to 



get Mfc and then integrating it using Eq. ( 19 1. The result 
is presented in Fig. [l] Frame (a) and (b) shows that the 
response, i.e., the transverse magnetization Mz which re- 
mains localized somewhere close to its initial value for all 
time. In other words, the response retains the memory 
of the breaking of the Z2 symmetry in transverse direc- 
tion by the polarized initial state through all later time, 
though the symmetry is respected by the driving over 
each complete cycle. The degree of symmetry-breaking 
is given by the long-time average of '■ 



lim 



1 



Mz{t)dt 



(23) 



Q is also a measure of non-adiabatic freezing - if a driv- 
ing were adiabatic, the resulting response would always 
follow the field (i.e., trace the instantaneous ground state 
value of the response) and thus would preserve the sym- 
metry of the Hamiltonian over a period. The maximum 
amplitude of oscillation of also determines the degree 
of freezing. 

It is clear from Fig. [TJthat for a given value of Fq, the non- 
adiabatic freezing Q is a strongly non-monotonic function 
of T. When the condition 



FoT 

h 



TT, 27r, Stt • 



(24) 



is satisfied the freezing attains a maximum (Q shows a 
peak) , as shown in Fig. [T](b) and [ijc) . Naively speaking, 
for a given Fq, if T is made larger, there is more time for 
system to react to the successive flips made, and hence 
the response is expected to be more adiabatic (smaller Q 
and bigger response amplitude). This classical intuition 
clearly does not hold in this case, as the freezing [Q and 
response amplitude) is strongly non-monotonic in T for 
a given Fg (Fig. Il]a). Strong maximal freezing of the 
entire many-body system {Q peaks) observed for isolated 
points in the parameter space is also a surprising non- 
classical feature of DMF, arising from coherent quantum 
dynamics [I]. 
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FIG. 1: The dynamical many-body freezing (DMF) behavior of the resulting response, (a) Variation of Mz with t for different 
p(p = ^) for To = 20. (b) Variation of Q with 2-r /T for To = 20. (c) Variation of Q with p for different Tq. Maximal 
Freezing are seen for integer p. (d) Magnetization after 100th and 1000th cycle at different To for T = 0.1. (ft = 1) 



In order to derive the extremal freezing condition ( 24 1 , 
we set T = and evaluate Mk{t — riT) as a function of 
n. Thus, we are basically looking at the start of every 
oscillation. Also, we assume that initially (at t — 0) the 
system was in the ground state for the transverse field at 
that moment. Thus, we set a = 1, /? = in Eq. (fl7|). 



use Eq. (21 1 there and obtain X-{nT) and a;+(nT)which 



is then substituted in Eq. (20). The result is 



where 



(25) 
(26) 



Mfc = Ak + Rk cos{2nujk + 4) 

Ak = cos^6'i + Qkfk, 
Rl = gl [fi + sm\2ei) sin^ ^1 sin^ cuk] , 

tanjfc = — sin(20i) sin^i sinwfc 

Jk 

with 

fk = sin(26'i) sin /^i coscj^. + sin(202) sm /i2 (27) 
gk = sin(2(?;>)sin(/i2)/(2sin^a;fc) = |C/i2|/(2sin^tjfc) 

and 



ujk = cos ^ [cos(/Lti + /X2) cos^ (p + cos(/Xi — 112) sin^ 



(28) 



From Eq. ( 20 ) and ( 25 1 we see, the non-adiabatic freezing 
parameter Q (Eq. 23 ) is given by 

g = -1 + ^ / Akdk (29) 



Now, for large Fq, from Eqs. (13 1 and (16) we get 
sin k 



vr 
■2 + 
FoT 

h 



Fn 



O 



cos k 



and 



(30) 



O 



The off-diagonal elements of the transfer matrix Ufe be- 
comes then, 



U 



12 



sin fi2 sin 20 



FoTX 2sinfc 



-0{1/Tl) 



(31) 



Hence, according to Eq. ( 16 ) if is an integral 



multiple of tt, Ufc becomes a Identity matrix up to 
terms I/Fq (since ^1 « —T^T/h for Fq ^ 1), and the 
system is found at the initial state (approximately) after 
each cycle. Note that the freezing occurs for any initial 
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state, irrespective of whether it is an eigenstate of the 
initial Hamihonian or not. It is also consistent with 
the Floquet picture of quasi-energy degeneracy (see e.g., 
[211 122| ) employed in explaining dynamical localization. 
According to Floquet theory, the above time evolution 
operator U^, which induces evolution from i = to 
t = T, should have the general form Ufc = e*'^'**, where 
Mk is a time independent Hermitian matrix (sharing 
same dimension and space as U^) |21| . with eigenvectors 
denoted by |/iifc) and \^2k) corresponding to eigenvalues 
Hik and /Z2fe, which are the Floquet quasi-energies. 
Now as we have shown, in our case Ufe tend to the 
Identity matrix up to term I/Fq in large Fq limit, which 
means both its eigenvalues e^^^'' , e^^^'' tend to unity 
for every k within the said approximation, resulting 
in a massive quasi-energetic degeneracy all over the 
many-body spectrum - the crux of DMF. Recently, 
another interesting manifestation of DMF is observed in 
periodically driven bosons in optical lattice with small 
amplitude driving across the tip of the Mott lobe [23] . 



The mechanism of DMF with rectangular driving can 
be visualized appealing to the simplicity of the driv- 
ing - it consists of dynamics driven by piecewise time- 
independent Hamiltonians. From Eq. (32) one can 



see, the dynamics (in the eigen-basis {K+Fq, /c)±)} of 
'Hfe(-)-Fo)) can be broken up into successive rotation of 
the basis by (due to successive flips of the transverse 
field) and intermediate accumulation of phases /ii_2 (due 
to intermediate waitings of length T/2). Clearly if one 
could adjust the intermediate phases /xi_2 such that their 
effect is nullified for all fc, in each cycle, then the system 
would return very closely to it's initial state after ev- 
ery cycle for any of the eigenstates {|(+Fo, A:)±)} as the 
initial state - the eigen-state of the initial Hamiltonian 
becomes Floquet states with degenerate quasi-energies 
(albeit with a difference in sign, which does not matter in 
this case). This happens, as explained in the paragraph 
following Eq. (32), when Fq is large and the condition 
for maximal freezing (Eq. 24) is satisfied. For small Fq, 
/ii 2 would retain strong fc-dependence, and hence this 
massive collapse of Floquet spectrum would not have 
been possible (see Fig. [l] d) . It is however worth not- 
ing that this simple picture of DMF cannot be extended 
in cases of continuously driven systems. For example, in 
the case of sinusoidal driving, the eigen-states of the ini- 
tial Hamiltonian do not tend to return to themselves as 
one approaches the DMF freezing peak - they retain a 
strongly fc-dependent period (fc is the quasi-momentum 
diagonalizing the initial Hamiltonian) of oscillation which 
actually diverges in the thermodynamic limit for certain 
modes as the DMF peak is approached. Freezing in that 
case is visualized as vanishing of the amplitude of oscil- 
lation of each fc mode, rather than "all fc-modes coming 
back to itself. It seems massive collapse of the entire 
DMF spectrum is a result of integrability of the model 
and simplicity of the driving. 



IV. LONG-LIVED SOLITARY OSCILLATION: 
THE SURVIVOR OF DESTRUCTIVE 
INTERFERENCES 

Analysis of Mz{t) shows that it is dominated by a 
distinct solitary oscillation in the long-time limit. The 
analysis of the response reveals sinusoidal oscillations 
of only two distinct timescales - one (denoted by ujq) 
matches with the driving period T (as expected), while 
the other, denoted by Tq (corresponding to frequency 
LJQ — 27r/TQ), depends on all the driving parameters. 
Tq can be much larger compared to T. In spite of the 
fact that the driving has a large amplitude and high fre- 
quency, and the system has several excitable energy lev- 
els, we observe only one distinct non-trivial frequency in 
the response. 

This can be understood as follows. The transverse 
magnetization Mz{t) (Eq. 19 1 at a time t is a super- 

For 



25) 



position of contributions for all fc's (Eq 
sufhciently large n, the argument {2nuj}^ + <5fe) in Eq.(25) 
will be large (so that its cosine will fluctuate very rapidly 
with fc) while Rk will remain relatively slowly varying. 
Thus the contributions from neighboring fc's will cancel 
out due to destructive interference (adding up with al- 
most same amplitude but rapidly varying phase) over any 
small intervals of fc, except those around the stationary 
points of Wfc (with respect to fc). In the neighborhood of 
its stationary points, Wk is expected to vary slowly with 
fc, and hence the contribution from different fc's within 
such neighborhood is expected to add up constructively. 
Elsewhere the contributions adds up destructively and 
can hence be ignored. Thus we may write 

/•f +e 

RkCos{2nujk+^k)dk Ri / cos{2niJk+^7T/2)dk 

(32) 

By Taylor expansion of uik about the stationary point, 
we can write. 



cos{2niiJk+STT/2) — cos{2nuj-^/2+ST^/2) cos I nC 









(nC 


\k--] 






r 21 





(33) 



where C = {d^uok/dk'^)k='K/2- This finally gives, 

M, (n) « Afo + cos(7ia;Q + S^/2) (34) 
V n 

where, 

ujQ = 2uj^/2 — 2cos^"^ {1 — cos^ (p[l — cos(^i -I- H2)]} 

(35) 

and 



Mo = -l 



2 r 



Ak dk, a = R 



"it/2 



Above arguments are quite generic, and variants of them 
can be found in different other contexts (see, e.g., [H]). 
Survival of few such distinct oscillations of very long 
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FIG. 2: Features of magnetization obtained by numerical integration of Eqs. (251 and (19 1, (a) Oscillations of two distinct 
time-scales are observed in Mz{t). In addition to the expected one with period T (matching to the driving), there is an 
additional longer time-scale prominently visible, (b) Fourier transform of M^{t) showing two timescales visible in (a): two 
distinct peaks at angular frequencies ujq and uiq are observed. The peak at ujq = 1.822 in the Fig, matches quite well with the 
estimation (1.815) from Eq. (351. While uiq « 2ti/T, where T = 0.1 is the driving period, (c) Long-time behavior of Mz{t) 
obtained numerically (points) and that given by Eq. fitting ( 34 1 (continuous curve) using is shown. The envelop corresponding 
to the l/\/n decay (Eq. 34 1 is visible, (d) Variation of Tq = 27r/ajQ with p for T — 0.1 obtained by Fourier transform. Tq 
tends to blow up (i.e., ujq vanishing up to order l/Fo) at integer values p - consistent with the observed maximal freezing at 
integer p. The points correspond to the values obtained from Fourier Transform and the continuous line is obtained from Eq. 
(351. (ft= 1) 



(compared to the driving period) time-scales was also ob- 
served in an infinite-range transverse Ising model driven 
periodically in time [53]. The results described in this 
section are manifestation of more general results re- 
garding periodically driven quantum many-body system 
([26]). 

We see from Eq. ( 35 ) , when freezing condition (p = in- 



teger) is satisfied, ujq vanishes for large Tq up to terms 
linear in I/Fq and a — )• 0, Mq — 1. Numerical calculation 
of Mz (using Eqs. (17 20)) is presented in Figj2j Discrete 
Fourier Transform of Mz(t) also shows two peaks corre- 
sponding to Wo and ujq. The value of luq obtained from 
there matches pretty well with the analytical expression 
in Eq.(l35l). 



Though our results are demonstrated for rectangular 
pulses, similar argument can be extended for other forms 
of periodic drivings. The only requirements for appear- 
ance of solitary oscillations (if they exist) are certain an- 
alytical properties of the response, continuity of the spec- 



trum, and long driving time. Hence such oscillations are 
expected to appear quite generically in many-periodically 
driven coherent many-body quantum system, but analyt- 
ical results might not be easy to extract in all cases. An 
extension of DMF for some other forms of periodic driv- 
ings may be achieved following [27] . 

The phenomena we discussed above are result of quan- 
tum coherence. Further investigations in this direction 
are likely to reveal many new phenomena (see. e.g. 
Ref.[28]). A natural open question is whether they are re- 
alizable in real experiments, in presence of the inevitable 
experimental imperfections existing within the present 
day setups. Such experimental realizations would also 
allow for exploring these phenomena in more generic non- 
integrable systems where accurate theoretical investiga- 
tions could be difficult. Experimental observation of the 
above phenomena might be possible within the frame- 
work of coherent quantum simulation using trapped ions 
and atom in optical lattice. In particular, DMF will 
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have clear signature even for very small systems con- 
sisting of few spins realized in the experimental systems 
above, since at the freezing peaks, all the momentum- 
modes freeze independent of system-size, whereas away 
from the peaks, considerable dynamics is expected for 
any system-size. Coherent simulation of transverse Ising 
Hamiltonian with time-dependent transverse field, which 
can be varied adiabatically, has also been realized exper- 
imentally. In layered linear Paul traps using ^'^^Yb^ ions 
[TT], and '^^Mg'^ ions [12], they realized transverse field 
Ising model where they can tune both spin-spin interac- 
tions and transverse field with time. 



V. SUMMARY 

We investigate dynamics of the transverse Ising chain 
under periodic instantaneous quenches of the transverse 
field. We make two interesting observations - 
(i) In the high amplitude (Fq 3> J) and fast quenching 
(T <C J) limits we observe Dynamical Many-body Freez- 
ing - we see that the driven system freezes close to its 
initial state, and the degree of freezing is a highly non- 
monotonic function of the pulse amplitude Fq and period 
T. The extremal freezing is observed for ToT/h = rnr 
(n = positive integers). At these freezing "peaks", the 
system remains frozen very strongly independent of its 
initial state. This freezing drastically contrasts the clas- 
sical notion of monotonic (with respect to the driving 
rate) freezing of a system under fast periodic driving - 
a faster driving would give it lesser time to react and 
hence would leave it more frozen. Quantum simulation 



of transverse Ising chain has already been realized ex- 
perimentally - the phenomenon should be amenable to 
experimental verification within the said set-up and sim- 
ilar others for quantum simulation. 

(ii) In the response dynamics, we observe emergence of 
a single, distinct timescale Tq (in addition to the time- 
scale of the driving) in the long-time limit. This dis- 
tinct oscillation decays much slower than other oscilla- 
tions, following a {n ^number of sweeps) envelop. 
Dominance of a single non-trivial frequency in the re- 
sponse is surprising, since the system is driven with pulses 
with high (compared to the intrinsic energy scale given 
by the spin-spin interaction J) amplitude and frequency. 
We show that this surviving time-scale represents oscil- 
lations of the non-local momentum modes lying within a 
neighborhood of a unique point in the momentum space 
(k = tt/2 here), where the contributions from the neigh- 
boring modes adds up constructively. For all other parts 
of the momentum space such interferences are destruc- 
tive, leading to mutual cancellation of oscillations of the 
neighboring modes. 
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